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We measure the ratio of cross sections, cr{pp Z + 2 b jets)/a{pp Z + 2 jets), for associated 
production of a Z boson with at least two jets with transverse momentum > 20 GeV and 
pseudorapidity < 2.5. This measurement uses data corresponding to an integrated luminosity 
of 9.7 fb“^ collected by the DO experiment in Run II of Fermilab’s Tevatron pp Collider at a center- 
of-mass energy of 1.96 TeV. The measured integrated ratio of 0.0236 ± 0.0032 (stat) ± 0.0035 (syst) 
is in agreement with predictions from next-to-leading-order perturbative QCD and the Monte Carlo 
event generators pythia and ALPGEN. 

PACS numbers: 12.38.Qk, 13.85.Qk, 14.65.Fy, 14.70.Hp 


Studies of Z boson production in association with a 
bottom and an anti-bottom quark provide important 
tests of the predictions of perturbative quantum chro¬ 
modynamics (pQCD) [H,[^. A good theoretical descrip¬ 
tion of this process is essential since it forms a major 
background for a variety of physics processes, including 
standard model (SM) Higgs boson production in associ¬ 
ation with a Z boson, ZH{H bb) Q, and searches for 
supersymmetric partners of the b quark [^. 


The ratio of Z -|- 6 jet to Z -|- jet production cross sec¬ 
tions, for events with at least one jet, has been previ¬ 
ously measured by the CDF [^, Q and DO 0-0 collabo¬ 
rations using Run II data. The ATLAS [l^ and CMS 111 
collaborations have also studied Z -|- 6 jet production at 
yi = 7 TeV. 
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This article presents the ratio of Z-l-2 b jets to Z-l-2 jets 
inclusive production cross sections and is an extension of 
the previous DO measurements utilizing similar event se¬ 
lections. The measurement of the ratio benefits from the 
cancellation of many systematic uncertainties, such as 
the uncertainty in luminosity and those related to lepton 
and jet identification, allowing a more precise compari¬ 
son with theory. The remaining systematic uncertainties 
arise from the differences between b jets and light jets. In 
the following, light-quark flavor (u, d, s) and gluon jets 
are referred to as “light jets”. The Z -b 2 6 jet produc¬ 
tion cross sections have been measured at CMS [l^ and 
ATLAS [10 at y/s = 7 TeV. The current measurement 
is based on the complete Run II data sample collected 
by the DO experiment at the Fermilab Tevatron pp col¬ 
lider at a center-of-mass energy of y/s = 1.96 TeV, and 
corresponds to an integrated luminosity of 9.7 fb“^. 


We first briefly describe the main components of the 
DO Run II detector IJ, 110 relevant to this analysis. The 
DO detector has a central trackin g sy stem consisting of a 
silicon microstrip tracker (SMT) 16| and a central fiber 
tracker (CFT), both located within a 1.9 T supercon¬ 
ducting solenoidal magnet, with designs optimized for 
tracking and vertexing at pseudorapidities |r7det| < 3 
and l^ydetl < 2.5, respectively [17|. A liquid argon and 
uranium calorimeter has a central section (CC) covering 
pseudorapidities |?7det| ^ 1-1, and two end calorimeters 
(EC) that extend coverage to hnst.l ~ 4.2, with all three 
housed in separate cryostats [l8| . An outer muon system. 
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at |77det| < 2, consists of a layer of tracking detectors and 
scintillation counters in front of 1.8 T toroids, followed 
by two similar layers after the toroids. Luminosity is 
measured using plastic scintillator arrays located in front 
of the EC cryostats. The trigger and data acquisition 
systems are designed to accommodate the high instanta¬ 
neous luminosities of Run II. 

This analysis relies on all components of the DO de¬ 
tector: tracking systems, the liquid-argon sampling 
calorimeter, muon system, and the ability to identify sec¬ 
ondary vertices 1^. The SMT allows for precise recon¬ 
struction of the primary pp interaction vertex and sec¬ 


ondary vertices [17|, [lOj. It also enables an accurate 


determination of the impact parameter, defined as the 
distance of closest approach of a track to the primary 
interaction vertex in the x-y plane. The impact param¬ 
eter measurements of tracks, along with reconstructed 
secondary vertices, are important inputs to the 6-jet tag¬ 
ging algorithm. 

Events containing Z bosons decaying to yy or ee 
are collected using triggers based on single electrons or 
muons. Eor the off-line selection requirements discussed 
below, the triggers have an efficiency of approximately 
100% for Z ^ ee and more than 78% for Z ^ yy decays 
depending on the transverse momentum of the muon. 
The Z -I- 2 jet sample requires the presence of at least two 
jets in the event, while the Z + 2 b jet sample requires 
at least two 6-jet candidates, selected using a 6-tagging 
algorithm (^ . 

An event is selected if it contains a pp interaction ver¬ 
tex, reconstructed from at least three tracks, located 
within 60 cm of the center of the DO detector along 
the beam axis. The selected events must also contain 
a Z boson candidate with a dilepton invariant mass 
70 < Mm < no GeV. 


Dielectron (ee) events are required to have two elec¬ 
trons of transverse momentum (pr) greater than 15 GeV 
identified through electromagnetic (EM) showers in the 
calorimeter. The showers must have more than 97% of 
their energy deposited in the EM calorimeter, be isolated 
from other energy depositions, and have transverse and 
longitudinal energy profiles consistent with that expected 
for electrons. At least one electron must be identified in 
the GG, with |?7det| < 1-1, and a second electron either 
in the GG or the EG, 1.5 < |??det| < 2.5. Electron can¬ 
didates in the CC are required to match central tracks 
or have a pattern of hits consistent with the passage of 
an electron through the central tracker. Electrons in the 
EGs are not required to have a track matched to them 
due to deteriorating tracking coverage for |?7det| > 2. Due 
to the lack of track requirement in EC regions we do not 
apply any opposite sign requirement for the dielectron 
events. 

The dimuon (yy) event selection requires two oppo¬ 
sitely charged muons detected in the muon system that 
are matched to reconstructed tracks in the central tracker 


withpT >15 GeV and |?7det| < 2. These muons must pass 
a combined tracking and calorimeter isolation require¬ 
ment discussed in detail in Ref. Q. Muons originating 
from cosmic rays are rejected by applying timing criteria 
using the hits in the scintillation counters and by limit¬ 
ing the measured displacement of the muon track with 
respect to the pp interaction vertex [2l| . 

A total of about 1.2 million Z boson candidate events 
are retained in the combined ee and yy channels with 
the above lepton selection criteria. The Z + 2 jet sample 
is then selected by requiring at least two jets in the event 
with > 20 GeV and I??-*®*! < 2.5. Jets are recon¬ 
structed from energy deposits in the calorimeter using 
an iterative midpoint cone algorithm with a cone 
of radius AR — ^/{A^p)^+{A^ = 0.5 where p is the 
azimuthal angle and y is the rapidity. Jet energy is cor¬ 
rected for detector response, the presence of noise and 
multiple pp interactions. We also correct the jet energy 
for the energy of those particles within the reconstruc¬ 
tion cone that is deposited in the calorimeter outside the 
cone (and vice versa) [i^ . 

To suppress background from top-antitop quark (it) 
production, events are rejected if the missing transverse 
energy is larger than 60 GeV, reducing the ti contamina¬ 
tion by a factor of two. These selection criteria retain an 
inclusive sample of 20,950 Z -|- 2 jet event candidates in 
the combined ee and yy channels. 

Processes such as diboson {WW, WZ, ZZ) produc¬ 
tion can contribute to the background when two lep¬ 
tons are reconstructed in the final state. Inclusive dibo¬ 
son production is simulated with the pythia [ 2 ^ Monte 
Carlo (MG) event generator. The Z jet, including 
heavy flavor jets, and tt events are modeled by ALP- 
GEN [ 2 ^, which generates hard sub-processes includ¬ 
ing higher order QCD tree level matrix elements, inter¬ 
faced with PYTHIA for parton showering and hadroniza- 
tion. The CTEQ6L1 [2^ parton distribution functions 
(PDFs) are used in all simulations. The cross sections of 
the simulated samples are then scaled to the correspond¬ 
ing higher-order theoretical calculations. For the diboson 
and Z -I- 2 jet processes, including the Z + bb signal pro¬ 
cess and Z + CC production, next-to-leading ord er ( NLO) 
cross section predictions are taken from MCFM [27|. The 
ti cross section is determined from NLOj-NNLL (next- 
to-next-leading log) calculations 28[. To improve the 
modeling of the px distribution of the Z boson, simu¬ 
lated Z-I-2 jet events are also reweighted to be consistent 
with the measured px spectrum of Z bosons observed in 
data [ 2 ^. 

These generated samples are processed through a de¬ 
tailed detector simulation based on GEANT [30|. To 
model the effects of detector noise and pile-up events, 
collider data from random beam crossings with the same 
instantaneous luminosity distribution as for data are su¬ 
perimposed on simulated events. These events are then 
reconstructed using the same algorithms as used for data. 
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Scale factors, determined from data using independent 
samples, are applied to account for differences in recon¬ 
struction efficiency between data and simulation. The 
energies of simulated jets are corrected, based on their 
flavor, to reproduce the resolution and energy scale ob¬ 
served in data (^ . 

The background contribution from multijet events, in 
which jets are misidentified as leptons, is evaluated from 
data. This is performed using a multijet-enriched sample 
of events that pass all selection criteria except for some of 
the lepton quality requirements. In the case of electrons, 
the multijet sample is obtained by inverting the shower 
shape requirement and relaxing other electron identifi¬ 
cation criteria, while for the muon channel, the multijet 
sample consists of events with muon candidates that fail 
the isolation requirements. The normalization of the mul- 
tijet background is determined from a simultaneous fit to 
the dilepton invariant mass distributions in different jet 
multiplicity bins. 

Figures [Hand [5] show the dilepton invariant mass and 
leading jet pr distributions in data compared to the ex¬ 
pectations from various processes. The dominant con¬ 
tribution comes from Z-|-light jet production. The non- 
Z + jet background fraction in the ee channel is about 
15%, and is dominated by multijet production. The 
muon channel has a higher purity with a background 
fraction of about 7%. 

This analysis employs a two-step procedure to deter¬ 
mine the 6-quark content of jets in the selected events. 
First, a 6-tagging algorithm is applied to jets to select 
a sample of Z -I- 2 jet events that is enriched in heavy 
flavor jets. After 6 tagging, the relative light, c, and 6- 
quark content is extracted by fitting templates built from 
a dedicated discriminant that provides an optimized sep¬ 
aration between the three components. 

Jets considered for 6-jet tagging are subject to a pre¬ 
selection requirement, called taggability, to decouple the 
intrinsic performance of the 6-jet tagging algorithm from 
effects related to the track reconstruction efficiency. For 
this purpose, the jet is required to have at least two asso¬ 
ciated tracks with pr > 0.5 GeV, the leading track must 
have pt > 1 GeV, and each track must have at least one 
SMT hit. This requirement has a typical efficiency of 
90% per jet. 

The 6-jet tagging algorithm is based on a multivariate 
analysis (MVA) technique (sij. This algorithm, MVAf,;, 
discriminates 6 jets from light-flavor jets utilizing the rel¬ 
atively long lifetime of the 6 hadrons when compared to 
their lighter counterparts Events with at least two 
jets tagged by this algorithm are considered. 

The MVAbi discriminant combines various properties 
of the jet and associated tracks to create a continuous 
output that tends towards unity for 6 jets and zero for 
light jets. Inputs include the number of secondary ver¬ 
tices and the charge track multiplicity, invariant mass of 
the secondary vertex (Msv), decay length and impact pa- 
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FIG. 1: (color online) The invariant mass in (a) Z —>■ 
pLpL and (b) Z ^ ee channels for data and background 
in events with a Z boson candidate and at least two 
jets before 6 tagging is applied. 
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rameter of secondary vertices, the multiplicity of charged 
tracks associated with them, and the Jet Lifetime Prob¬ 
ability (JLIP), which is the probability that tracks as¬ 
sociated with the jet originate from the pp interaction 
vertex Events are retained for further analysis if 

they contain at least two jets with an MVAf,/ output 
greater than 0.15. After these requirements, 241 Z-\-2 jet 
events are selected with at least two 6-tagged jets, where 
only the two highest px tagged jets are examined in the 
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FIG. 2: (color online) The leading jet px in the (a) 

Z —^ pp and (b) Z ^ ee channels for data and back¬ 
ground in events with a Z boson candidate and at least 
two jets before b tagging is applied. 
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FIG. 3: (color online) The one dimensional projection 
onto (a) the highest-pT jet and (b) the second highest- 
Pt jet Hmjl axis of the 2D fit. The distributions of 
the b, c, and light jets are normalized by the fractions 
found from the fit. 


analysis and the electron and muon channels are com¬ 
bined. The efficiency for tagging two b jets in data is 
33%. In the MC correction factors are applied to account 
for differences with data [i^. The background contami¬ 
nation from diboson, multijet, and top production after 
6-tagging, for the electron and muon channels combined 
are 8%, 2% and 15% respectively. 

To determine the fraction of events with 2 6 jets, a 
dedicated discriminant, Dmjl, is employed [ 1 , | 32 |. It 
is a combination of the two most discriminating MVA;,; 


inputs, Msv and JLIP: Dmjl = 0.5 x (Msv/(5 GeV) — 
ln(JLIP)/20). 

To measure the fraction of events with different jet 
flavors in the selected sample, we count the number of 
events as a function of the Dmjl of the two leading jets 
A^(T*mjLi) T^mjl 2 ) and then perform a two dimensional 
binned maximum likelihood fit to that distribution. The 
data sample with two heavy-flavor-tagged jets is fitted to 
templates consisting mainly of 2 6-jet, 2 c-jet, and light 
flavor jet events, as obtained from ALPGEN-|-pythia sim- 
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ulated samples. We also compared the shapes of the tem¬ 
plates from SHERPA simulated samples and found the 
templates to be consistent for the two models. Before the 
fit, all non-Z + jet background contributions, estimated 
from simulated samples after the MVA;,; requirement, are 
subtracted from the data leaving 180 Z + 2 jet events in 
the combined ee and fifj, channel. Next, we measure the 
jet-flavor fractions in the dielectron and dimuon sam¬ 
ples combined, yielding the 2 b jet flavor fraction (/{,&) 
of 0.64 ± 0.08 (stat.) and the 2 c jet flavor fraction of 
0.32 ± 0.08 (stat.). Figure |3] shows the one dimensional 
projection onto the highest-pr jet and the second highest- 
Pt jet Dmjl axis of the 2D fit. 

The fraction of 2 6 jets measured in the heavy fla¬ 
vor enriched sample is combined with the corresponding 
event acceptances to determine the ratio, R, of the cross 
sections. 


^ a{Z + 2b jets) Ntt fbb Anci 
a{Z + 2 jets) Nind etag ^ Ab 

where A^inci is the total number of Z+ 2 jet events before 
the tagging requirements, Nf,/, is the number of Z -|- 2 jet 
events used in the Dmjl fit, fbb is the extracted 2 b 
jet fraction, and is the overall selection efficiency 
of Dmjl for 2 b jets that combines the efficiencies for 
taggability and MVA;,; discriminant. Both Aind and Nbb 
correspond to the number of events that remain after 
the contributions from non-Z jets processes have been 
subtracted from the data. 

The pseudorapidity acceptance for electrons and 
muons is different. In order to quote a combined ra¬ 
tio for the two channels, we correct to a common lepton 
acceptance as follows. The detector acceptances for the 
inclusive jet sample and 2 b jets are determined from MC 
simulations in the kinematic region that satisfies the px 
and r] requirements for leptons and jets. For the Abb and 
Ainci calculations, we apply selections for both the elec¬ 
tron and muon channels for the fiducial region for the 
events with two jets and two leptons defined as: 

> 20 GeV and < 2.5, 

p^ > 15 GeV and \/\ < 2. (2) 

The resulting ratio of the two acceptances is measured 
to be Amd/Abb = 1-09 ± 0.02 (stat). 

Using Eq. o, we obtain the ratio of the Z -|- 2 6 jet 
cross section to the inclusive Z -|- 2 jet cross section in the 
combined pp and ee channel to be 0.0236±0.0032 (stat). 

Several systematic uncertainties cancel when the ra¬ 
tio ct(Z ±26 jets)/cr(Z ± 2 jets) is measured. These in¬ 
clude uncertainties on the luminosity measurement, lep¬ 
ton trigger efficiency, and lepton and jet reconstruction 
efficiencies. The remaining uncertainties are estimated 
separately for the integrated result. The largest system¬ 
atic uncertainty of 13.7% comes from the uncertainty on 


the shape of the Dmjl templates used in the fit includ¬ 
ing that due to MG statistics of the samples used to con¬ 
struct the templates. The shape of the templates may be 
affected by the choice of the 6 quark fragmentation func¬ 
tion 33| , the background estimation, the difference in the 
shape of the light jet MG template and a template de¬ 
rived from a light jet enriched dijet data sample, and the 
composition of the charm states used to determine the 
charm template shape Q. It also includes uncertainties 
on production rates of different hadrons and uncertain¬ 
ties on branching fractions. These effects are evaluated 
by varying the central values by the corresponding un¬ 
certainties, one at a time. The entire analysis chain is 
checked for possible biases using a MC closure test and 
no significant deviations are observed. The next largest 
systematic uncertainty of 5.5% is due to the 6-jet iden¬ 
tification efficiency. The uncertainty on 6-jet energy cal¬ 
ibration is 2.6%; it comprises the uncertainties on the 
jet energy resolution and the jet energy scale. For the 
integrated ratio measurement, these uncertainties, when 
summed in quadrature, result in a total systematic un¬ 
certainty of 14.9%. For the integrated ratio we obtain 


R = 0.0236 ± 0.0032 (stat) ± 0.0035 (syst). 


(3) 


To check the stability of the result, the ratio is remea¬ 
sured using a looser (tighter) MVAf,; selection with the 
lower limit on the MVAf,; output of > 0.10(> 0.225). The 
looser selection provides increased data statistics and the 
tighter selection yields a 2 6 enriched sample. The new 
and default ratios are found to be in agreement within 
uncertainties of about 4%. 

To validate the tt background estimation, we reduce 
the contribution of tt events by rejecting events where the 
scalar sum of all jet px values is more than 130 GeV. This 
selection reduces the tt fraction by an additional factor 
of 2 with a signal efficiency of 80%. The new and default 
ratios are found to be in agreement within systematic 
uncertainties. 

In Table m we present the ratio of integrated cross sec¬ 
tions, <t{pp — >• Z ± 2 6 jet)/CT(pp —>■ Z ± 2 jet), in the 
fiducial region defined in Eq. ©■ The ratio is compared 
to predictions from NLO QCD calculations and two MC 
generators, pythia and ALPGEN. The NLO predictions 
use the MSTW2008 PDF set 341 using mcfm with cen¬ 
tral values of renormalization and fragmentation scales 
Pr = Pf = Mz- Uncertainties are estimated by varying 
Pr and p f together by a factor of two, and are about 15%. 
ALPGEN generates multi-parton final states using tree- 
level matrix elements. When interfaced with pythia, 
it employs an MLM scheme 3^ to match matrix ele¬ 
ment partons with those after showering in pythia, re¬ 
sulting in an improvement over leading-logarithmic ac¬ 
curacy. The measured ratio is in reasonable agreement 
with MCFM NLO calculations considering the uncertain¬ 
ties on the data and theory. 








TABLE I: The ratio of integrated cross sections, a{pp Z + 2 b jet)/a{pp Z + 2 jet) together with statistical 
uncertainties (5stat) and total systematic uncertainties ((5syst)- The column (5tot shows the total experimental 
uncertainty obtained by adding Jgtat and Jgyst in quadrature. The last three columns show theoretical predictions 
obtained using NLO QCD with scale uncertainties and two MC event generators, pythia and ALPGEN. 


cr{pp Z + 2 b jet)/(T(pp Z + 2 jet) 


DSjtSj lll(5stat i (^syst *5tot 

NLO QCd(mSTW) pythia ALPGEN 

(2.36 ±0.32 ±0.35) x 10"^ 0.47x10"^ 

(1.76 ± 0.26) X 10“^ 2.42x10"^ 2.21x10"^ 


In summary, we report the measurement at the Teva- 
tron of the ratio of integrated cross sections, a{pp —>■ 
Z + 2 b jet)/(T(pp Z + 2 jet), for events with Z ^ U 
in a restricted phase space of leptons with > 15 GeV, 
|? 7 ^| < 2.0 and with two jets limited to > 20 GeV 

and < 2.5. Measurements are based on the full 

data sample collected by the DO experiment in Run II 
of the Tevatron, corresponding to an integrated luminos¬ 
ity of 9.7 fb~^ at a center-of-mass energy of 1.96 TeV. 
The measured integrated ratio of 0.0236 ±0.0032 (stat)± 
0.0035 (syst) is in agreement with the theoretical predic¬ 
tions within uncertainties. 

We thank John Gampbell and Doreen Wackeroth for 
valuable discussions, and the staffs at Fermilab and col¬ 
laborating institutions, and acknowledge support from 
the Department of Energy and National Science Foun¬ 
dation (United States of America); Alternative Ener¬ 
gies and Atomic Energy Commission and National Cen¬ 
ter for Scientific Research/National Institute of Nuclear 
and Particle Physics (France); Ministry of Education 
and Science of the Russian Federation, National Re¬ 
search Center “Kurchatov Institute” of the Russian Fed¬ 
eration, and Russian Foundation for Basic Research 
(Russia); National Council for the Development of Sci¬ 
ence and Technology and Carlos Chagas Filho Foun¬ 
dation for the Support of Research in the State of 
Rio de Janeiro (Brazil); Department of Atomic En¬ 
ergy and Department of Science and Technology (In¬ 
dia); Administrative Department of Science, Technol¬ 
ogy and Innovation (Colombia); National Council of Sci¬ 
ence and Technology (Mexico); National Research Foun¬ 
dation of Korea (Korea); Foundation for Fundamen¬ 
tal Research on Matter (The Netherlands); Science and 
Technology Facilities Council and The Royal Society 
(United Kingdom); Ministry of Education, Youth and 
Sports (Czech Republic); Bundesministerium fiir Bildung 
und Forschung (Federal Ministry of Education and Re¬ 
search) and Deutsche Forschungsgemeinschaft (German 
Research Foundation) (Germany); Science Foundation 
Ireland (Ireland); Swedish Research Council (Sweden); 
China Academy of Sciences and National Natural Science 
Foundation of China (China); and Ministry of Education 
and Science of Ukraine (Ukraine). 


[1] J. M. Campbell, R. K. Ellis, F. Maltoni, and S. Willen- 
brock, Phys. Rev. D 73, 054007 (2006). 

[2] F. F. Cordero, L. Reina, and D. Wackeroth, Phys. Rev. 
D 78 , 074014 (2008). 

[3] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
109 , 121803 (2012); T. Aaltonen et al. (CDF Collabora¬ 
tion), Phys. Rev. Lett. 109 , 111803 (2012). 

[4] V. M. Abazov et al. (DO Collaboration), Phys. Lett. B 
693 , 95 (2010); T. Aaltonen et al. (CDF Collaboration), 
Phys. Rev. Lett. 105 , 081802 (2010). 

[5] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
79 , 052008 (2009). 

[6] A. Abnlencia et al. (CDF Collaboration), Phys. Rev. D 
74, 032008 (2006). 

[7] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
94, 161801 (2005). 

[8] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D 
83 , 031105 (2011). 

[9] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D 
87 , 092010 (2013). 

[10] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 706 , 
295 (2012). 

[11] S. Chatrchyan et al. (CMS Collaboration), JHEP 06 
(2012) 126. 

[12] S. Chatrchyan et al. (CMS Collaboration), JHEP 06 
(2014) 120. 

[13] G. Aad et al. (ATLAS Collaboration), JHEP 10 (2014) 
141. 

[14] V.M. Abazov et al. (DO Collaboration), Nncl. Instrnm. 
Methods Phys. Res. A 565 , 463 (2006). 

[15] M. Abolins et al., Nncl. Instrnm. Methods Phys. Res. A 
584 , 75 (2008). 

[16] R. Angstadt et al., Nncl. Instrum. Methods in Phys. Res. 
A 622 , 298 (2010); S. Ahmed et al, Nucl. Instrnm. Meth¬ 
ods Phys. Res. A 634 , 8 (2011). 

[17] We nse a standard right-handed coordinate system. The 
nominal collision point is the center of the detector with 
coordinates (0, 0, 0). The direction of the proton beam is 
the positive +z axis. The +x axis is horizontal, point¬ 
ing away from the center of the Tevatron ring. The +y 
axis points vertically upwards. The polar angle, 9, is de¬ 
fined snch that 0 = 0 is the +z direction. The rapidity 
is defined as y = —ln[{E + pz)/{E—pz)], where E is 
the energy and pz is the momentnm component along 
the proton beam direction. Pseudorapidity is defined as 
p = — ln(tan|). (p is defined as the azimnthal angle in 



9 


the plane transverse to the proton beam direction. Also, 
?7det and (pdet are the pseudorapidity and the azimuthal 
angle measured with respect to the center of the detector. 

[18] S. Abachi et al. (DO Collaboration), Nucl. Instrum. Meth¬ 
ods Phys. Res. A 338, 185 (1994). 

[19] The primary pp interaction vertex is that found to be 
the most likely collision point, among possibly several 
collisions within a specific beam crossing, from which se¬ 
lected objects emanate. The algorithm for defining it can 
be found in [^ . 

[20] V. M. Abazov et al. (DO Collaboration), Nucl. Instrum. 
Methods Phys. Res. Sect. A 763, 290 (2014). 

[21] V. M. Abazov et al. (DO Collaboration), Nucl. lustrum. 
Methods Phys. Res. Sect. A 737, 281 (2014). 

[22] G. C. Blazey et al., arXiv:hep-ex/0005012 

[23] V. M. Abazov et al. (DO Collaboration), Nucl. Instrum. 
Methods Phys. Res. Sect. A 763, 442 (2014). 

[24] T. Sjostrand, S. Mrenna, and P. Skands, J. High Energy 
Phys. 05 (2006) 026. Version 6.409 with Tune A was used. 

[25] M. L. Mangano et al, J. High Energy Phys. 07 (2003) 


001. Version 2.11 was used. 

[26] J. Pumplin et al., J. High Energy Phys. 07 (2002) 012. 

[27] J. M. Campbell and R. K. Ellis, Phys. Rev. D 60, 113006 
(1999); %bid. 62, 114012 (2000); ibid. 65, 113007 (2002). 

[28] U. Langenfeld, S. Moch, and P. Uwer, Phys. Rev. D 80, 
054009 (2009). 

[29] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 
100 , 102002 (2008). 

[30] R. Brun and F. Carminati, CERN Program Library Long 
Writeup W5013 (1993). 

[31] A. Hoecker et al., TMVA: Toolkit for Multivariate Data 
Analysis. PoS, ACAT:040, 2007. 

[32] V. M. Abazov et al. (DO Collaboration), Phys. Lett. B 
718, 1314 (2013). 

[33] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D 
84, 032004 (2011). 

[34] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, 
Eur. Phys. J. C 63, 189 (2009). 

[35] F. Caravaglios et al., Nucl. Phys. B 539, 215 (1999). 


